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Abstract 

We consider a phenomenologically viable SO(IO) grand unification model which 
allows perturbative calculations up to the Planck scale or the string scale. We use 
a set of Higgs superfields 10 + 16 + 16 + 45. In this framework, the data fitting of 
the charged fermion mass matrices is re-examined. This model can indeed reproduce 
the low-energy experimental data relating the charged fermion masses and mixings. 
As for the neutrino sector, we take the neutrino oscillation data as input data to 
construct right-handed Majorana neutrino mass matrix and get a prediction for the 
physics related to the right-handed neutrinos, e. g. the leptogenesis and for the proton 
decay. We propose two kinds of phenomenologically viable model, quoted as Model 
1 and Model 2. We show that one of the models (Model 2) is consistent with all 
experimental constraints. 
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1 Introduction 



The supersymmetric (SUSY) grand unified theory (GUT) provides an attractive implica- 
tion for the understandings of the low-energy physics. In fact, for instance, the anomaly 
cancellation between the several matter muhiplets is automatic in the GUT based on a 
simple gauge group, since the matter muhiplets are unified into a few muhiplets, the ex- 
perimental data supports the fact of unification of three gauge couplings at the GUT scale 
Mq ~ 2 X 10^^ [GeV] assuming the particle contents of the minimal supersymmetric Stan- 
dard model (MSSM) and also the right-handed neutrino appeared naturally in the SO(IO) 
GUT provides a natural explanation of the smallness of the neutrino masses through the see- 
saw mechanism [2]. While the SUSY GUT has several nice feature described above, it should 
be available up to a scale, several order of magnitude above the GUT scale in order to open 
up a window for really a GUT picture. In other words, we require to have perturbative de- 
scription of the SUSY GUT up to a cutoff scale, say the Planck scale or the string scale. Once 
we adopt to have such a landscape, we can not include a large dimensional representations 
like 126 + 126 of SO(IO) as a Higgs since it causes a strongly asymptotic non-free behavior. 
In addition, a free-field heterotic string theory can never give a 126 + 126 representations 
[3]. Thus we consider a model based on a set of Higgs with small dimensional representa- 
tions and we explore the minimal choice of the Higgs superfields which can accommodate 
the low-energy experimental data. As for the neutrino sector, we take input the neutrino 
oscillation data and get a prediction for the physics related to the right-handed neutrinos, 
e. g. the leptogenesis, etc. 

Here it needs some comments for the alternative approaches are in order. There could be 
another possibility which keep the coupling constant perturbative up to some cut-off scale 
(that could be the string scale or the Planck scale) by including the threshold corrections as 
shown in Refs. [1]. Those approaches are also interesting as the masses of the colored Hig- 
gsinos can be raised enough to avoid too rapid proton decay without fine-tuning. However, 
at least, in this paper, we keep the simple unification picture within the framework of the 
MSSM particle content at low energy. 

Now we think about SO(IO) problems again. A very concrete problem we can ask and 
solve is the foUowings. Is there an SO(IO) model that has: 

(1) Gauge couplings unification similar to MSSM RG analysis. 

(2) Perturbative unification up to the string or the Planck scale. 

(3) Fit all known phenomenology. 

We think this is a very constrained requirement already. To demand perturbative uni- 
fication up to the string or the Planck scale will require the model is pretty "minimal" for 
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T(R) 


10 


1 


16 
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45 
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54 


12 


120 


28 


126 


35 


210 


56 



Table 1: List of the Dynkin index for the SO(IO) irreducible representations up to the 210 
dimensional one 



sure. There are two lines of minimal SO(IO) model. One uses higher dimensional represen- 
tation like 126 and renormalizable [5] [6] [7j [8], and the other does not use representations 
higher than 54 representation but non- renormalizable operators [9]. In order for satisfying 
the requirement (2) in the above list of questions, we investigate the allowed sets of Higgs 
superfields by considering their contributions to the beta function coefficient. In general, we 
have the RG equation for the unified gauge coupling ac, 

1 _ 1 b 

When we set the condition that the coupling constant «g blows up at a scale /i = A we have 
the following condition, 

1 27r , / A \ 

^— = 0, i.e. , = Mog — . 2 

«G (A) «G (Mg) v Mg/ 

If we want to have a perturbative description up to the (reduced) Planck scale Mpi = 2.4 x 10^^ 
GeV, we get the upper bound on the coefficient b as 

b < 34. (3) 

In general, the coefficient b can be written as 

b = J2 T(R)-3x ^ . (4) 

chiral multiplet vector multiplet 

When we extract the 3 generation matter contributions we have 
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b = 2^ T(R) + 3x2-3x8. 

Higgs multiplet 



(5) 



So the maximal value for the sum of T(R) is given by 



T(R) < 34 - 6 + 24 = 52. (6) 



Higgs multiplet 



For example, a set of Higgs multiplets {lO © 10' © 16 © 16 © 45} is allowable from pertur- 
bative argument, 



Y T{R) = 1 x2+ 2 x2+ 8 = 14 < 52. (7) 



Higgs multiplet 10 010' 16ffil6 45 



In the former minimal SO(IO) model including Higgs superpotential [10], the set of Higgs 
multiplets is given by 10 © 126 © 126 © 210, thus, in this case, we have a very huge size of 
the beta function coefiicient, b = 109. In this case, the cutoff scale obtained is very close 
to the GUT scale 



So this model is severe but not excluded since we have nothing definite beyond Mg and 
even Mg itself. Its highly predictivity is deserved to study further. However, in this letter, 
we consider that the existence of the string or the Planck scale is an reality that one has 
to accept. So it is actually necessary to include that. It is much more important to see if 
one can get perturbative unification up to the string or the Planck scale even when higher 
dimensional operators are included. While there are a priori many choices for the Higgs 
superfields, the predictivity and the perturbativity may pick up some representations. And 
we explore the minimal choice of the Higgs superfields which can accommodate the low- 
energy experimental data. As for the neutrino sector, we take input the neutrino oscillation 
data and get a prediction for the mass matrix of the right-handed neutrinos. 

2 Fermion mass matrices 

We can consider some possible choices for Higgs representations to be introduced in our 
SO(IO) model. In any case, they should accomplish the following tasks: (1) it can break 
the SO(IO) gauge symmetry down to the Standard model one. (2) it can reproduce all the 
fermion mass matrices being realistic. For the task (1), we give some examples, 45 + 54, 
16 + 16 + 45 etc. Considering the second task and the fact that top Yukawa coupling is 
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of order one, it is necessary to introduce, at least, one 10 Higgs representation which can 
provide a renormalizable Yukawa coupling. Note that 126 Higgs representation is excluded 
by our criteria of the perturbative unification, since 126 + 126 contributes &Higgs = 70. 
Furthermore, in order to incorporate right-handed neutrino Majorana masses, we need 16 
Higgs, and consider the superpotential 

W = j^Y^16d6,16Hl6H, (9) 

where M is the cutoff scale of our model, the Planck scale or the string scale, for example. 
In this paper, we consistently assume that only the singlet in 16h (t^h) and 16h (j^h) can 
achieve their VEVs and the doublets never get their VEVs. This assumption is essential 
to write down the GUT mass relations for the charged fermions. Then we may define our 
"minimal model" as the one which satisfies all the above requirements and contributes &Higgs 
as small as possible. The most reasonable candidate is the choice with 10 + 16 + 16 + 45 
Higgs multiplets. 

With these Higgs multiplets, the superpotential relevant to the fermion mass matrices is 
given by (up to dimension 5 terms) 

W = Yl^WaGA^H + ^^^16,16,10^^45^ (10) 

+ (7,16,1645^ + /1^1646 IOh (11) 
+ gl61645H + 16 16 10 + /161610// (12) 
+ MiolO^ + Mi6l6//T6// + M4545I + X16h1Qh45h , (13) 

where the Yukawa coupling matrices Yio is symmetric, while 145 is antisymmetric. Here, an 
extra vector like {16+ 16}-multiplet with no subscript represents an extra matter multiplet, 
whose matter parity is assigned to be odd, so that the Yukawa coupling given in Eqs. ffTOl) - 
( fT3|) becomes invariant under the matter parity. 

Note that the 16j and 16 do mix with each other. So, in general, the Yukawa couplings 
can be written as 

_ r,,10// {h,/2)10H (9^/2)45 

W = ( 16„ 16, 16)1 ihj/2)10H hlOn (^/2)45 

{gj/2)45H {~g/2)45H flOn 

_ 0{Mw) 0{Mw) 0{Mg) 
( 16,, 16, 16 ) I 0{Mw) 0{Mw) 0{Mo) I I 16 1 , (14) 



0(Mg) 0(Mg) 0(M, 
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where we provided the VEVs to the Higgs, (10//) ~ C(Mvk) and (45//) ~ 0{Mq) and 
this mass matrix has only one weak scale mass (9(Mw) and the remaining two are of order 
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the GUT scale 0{Mq). Then three lighter modes in the multiplet 16q, = (16^, 16, 16)^ 
(a = 1, 2, 3, 4, 5) are identified with the 16-multiplet including the usual quarks and leptons. 
The hght modes (169 — l!2,3) in the mass matrix of Eq. (fT4l) . which contains the 
usual quarks and leptons can be written as the foUowing linear combinations (up to some 
normalization factor): 

16; = -gl6i + gil6, 
I62 = -5(3162 + 5-1163 , 

16'3 = -52I63 + 51I62. (15) 

The remaining two eigenstates become heavy of order Mq, so the low energy effective theory 
can indeed become the MSSM without exotics. Then the extra vector-like matter © g^, 
ii^O) i'i, etc. would decouple from the low-energy effectve field theory. After giving a VEV 
to the 45 H Higgs multiplet, we have the dimension 5 neutrino mass operator Eq. iQ and 
also the dimension 6 operator contributing to the charged fermion mass matrices 

W=^ Y^16dQjlQHlQH + (1^) 1646,10^45^. (16) 

In the followings, we consider two cases for the description of the charged fermion mass 
matrices. One is to consider the situation in which the second term in Eq. ffTOl) has larger 
contribution than Eq. f|T6|) . Another one is the opposite case in which the dominant term is 
Eq. ffTBl) rather than the second term in Eq. flTUl) . 

In the former case, Eq. flTU]) corresponds to the Dirac mass matrices of quarks and leptons. 
Note that this example is obviously unrealistic since it predicts the Kobayashi-Maskawa 
matrix being unity. This is because the Higgs doublets in Eq. (ITO!) is the same, and as result 
the up-type quark mass matrix is proportional to the down-type quark mass matrix. One 
way to avoid this problem is to introduce new 10 Higgs and a symmetry which allow the 
superpotential such as 

W = ^1^16,16,101 + 16,16, 10245^, (17) 

where lOi and IO2 are two Higgs multiplets of 10 representation. Since 10 x 45 D 
10 + 120 + 320, this system is effectively the same as the one with 10 + 120 Higgs multi- 
plets. 

In the latter case, the system is completely equivalent to the minimal SO(IO) model 
which uses 10 + 126 Higgses w. r. t. the charged fermion sector since 10 X 452 D 126 + 126. 
In this case we need only one 10 multiplet to reproduce realistic fermion mass matrices. 
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Model 1 



In the foUowing notation, we use 120 Higgs, for simplicity. This model has been studied in 
jllj . in detail, and we give a very brief review of this model. 

The Yukawa couplings relevant to the Dirac mass matrices are given by 

W = YiimilG.lOH + 17^1646,120^^, (18) 

where Yiq and Y120 are symmetric and anti-symmetric, respectively. Number of free parame- 
ters in the Yukawa matrices is found to be 3 + 3 x 2 = 9 in total. Both of the Higgs multiplets 
lOj^ and 12Qh include a pair of Higgs doublets in the MSSM decomposition. At low-energy 
after the GUT symmetry breaking, the superpotential leads to 

+ iYllH-,, - Wll,H-,,)N,l, + {YllHf, - 3r/^o^i2o)e-^., (19) 

where Hiq and Hi2q correspond to the Higgs doublets in 10 an 120 Higgses, which originate 
lOi and IO2 in the original superpotential Eq. (fT7|) . The factor 3 in the lepton sector is the 
results from the VEV of 45 Higgs in the B — L direction. Providing VEVs for all the Higgs 
doublets, the Dirac mass matrices are obtained as 

W^ass = Mi^^u'rg, + M'Jd'rg, + M'^N.ij + Ml^e'ri,, (20) 

where 

Mu = cioMio + C120M120, 
Ma = M10 + M120, 

Md = CioMio - 3Ci2oMi20, 

Me = M10-3M120. (21) 

Here Cio = (-f^io) / {Hfo) and C120 = {H120) I {-^120) complex parameters in general. Now 
number of free parameters in terms of mass matrices are found tobe3 + 3x2 + 2x2 = 13. 
This relation leads to a relation among the mass matrices of up- and down-type quarks and 
charged lepton such as 

Me = Cd {Md + kMu) , (22) 

^ In general, there could be another contribution to the mass matrices from bi-doublet (1,2,2) C 120 
under G422 = SU(4) x SU(2)i x SU(2)ij. In this paper, we assumed that it does not contribute to the fermion 
mass matrices, just for simphcity. 
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where 



3cio + C120 

Cd = , 

ClO — C120 

«: = -7r4 • (23) 

OCio + C120 

Note that this GUT relation holds at the GUT scale. By the same methods in [7], we will 
solve the relation and find mass matrices compatible to the low-energy data of the fermion 
mass matrices. 

Without loss of generahty, we can begin with the basis where is real and diagonal, 
Mu = D u- Here we assume to be the hermitian matrix, and the number of parameters 
are reduced to 3 + 3 + 4 = 10. Therefore it can be diagonalized by an unitary matrix, = 
VKuDaV^y^. Considering the basis-independent quantities, tr(Me), tr (Mf ) and det (Mg), 
and eliminating Cd, we obtain two independent equations, 



tr(Me) \ tr(Me') 



me + + nir J ml + + 

tr(Me) \^ det(M;) 



(24) 
(25) 



where Mg = V^MDdV^y^ + kD^- With input data of six quark masses, three angles and 
one CP-phase in the CKM matrix and three charged lepton masses, we can solve the above 
equations and determine k and Cd- The original basic mass matrices, Mio and M120, are 
described by 

Mio = ^^V^M,D,Vl^+''-^ (26) 

M120 = ^-^VKmDdVly,-"-^ Du. (27) 

Now Mio and Mi2n are completely determined with the solutions Cd and k determined by 
the GUT relation. □ It should be checked that Y]*2o should be much smaller than one, since 
^120 — I ^ ■ the next section, we will perform a numerical analysis by using the 

same method as 



^ Although the solution is found in the GUT relation are appUed at the electroweak scale. The 
analysis incorporating renormalization effects was done as in [TJ . 
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Model 2 



The Yukawa couplings relevant to the Dirac mass matrices are given by 



W = F/o^ieae.lOH + F/2'6l6il6,126H, (28) 

where both Yio and Y126 are symmetric, and the number of free parameters in the Yukawa 
matrices is found to be 3 + 6 x 2 = 15 in total. In this model, the GUT relation is the 
same as the one in Eq. (|22l) . But in this case, the coefficients cio and Ci26 are complex 
numbers, and then the analysis for finding a solution is somewhat restrictive. Notice that 
the right-handed neutrino Majorana mass matrix is completely free in this model as depicted 
in Eq. ([9]), thus we can reproduce the neutrino oscillation data. But even in this case, since 
we know the neutrino Dirac mass matrix, we can get a prediction for the quantity relating 
the leptogenesis. 

Now as in the same way in Model 1, considering the basis- independent quantities, tr[MjMe], 
tr[(MtMe)2] and det[M]Me], and eliminating we obtain two independent equations. 



tr[M/Me] V _ tr[(MeWe)'] 



m? + + ml j mj + mf, + 



tr[MeVe] A det[MeVe 



+ m'j, + m'i j ml mj, mi 



(29) 



(30) 



where Me = U* W + nD^- With input data of six quark masses, three angles and one 
CP-phase in the CKM matrix and three charged lepton masses, we can solve the above 
equations and determine k and but five parameters, the phase of q [7] and the phases 
in U, except for the CKM phase are undetermined. The original basic mass matrices, Mio 
and Mi26, are described by 

M,o = ^^^^V^^D,V^^ + \^^D^, (31) 

By giving a solution for |q| and k, we can determine all the Dirac mass matrix as a function 
of a phase a. As a result, we can obtain a prediction for the leptogenesis parameters. 



''in this paper, we ignore the additional phases in U except for the KM phase S, and thus hereafter we 
denote U just as Vkm- 
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3 Numerical analysis and results 



Now we are ready to perform all the numerical analysis. We first sliow the solution in 
Model 1. In the following analysis, we assume all the charged fermion mass matrices to be 
hermitian. We take input the absolute values of the fermion masses at Mz as foUows (in 
GeV) [E]: 

mu = 0.00233, = 0.677, m* = 176, 
rud = 0.00469, m, = 0.0934, rrib = 3.00, 

TTie = 0.000487, m/, = 0.103, = 1.76. (33) 

Here the signs of the input fermion masses have taken as {mu,mc,mt) = (—,—,+) and 
{md,ms,mb) = (—,—,+). And for the CKM mixing angles and a CP-violating phase in 
the " Standard" parameterization, we input the center values measured by experiments as 
foUows: 



si2 = 0.2229, S23 = 0.0412, S13 = 0.0036, S = n/3 [rad]. (34) 

Since it is very difficult to search all the possible parameter region systematically, we present 
only reasonable results we found. In the following, we show our analysis in detail in the case 
tan /3 = 50. After the RG running, we get the values at the CUT scale are use them as input 
parameters in order to solve the GUT relation of Eq. (122|) . 

Now for simplicity, we assume the matrices Mio and M120 in Eq. (12T|) to be hermitie. 
Thus the coefiicients Cd and k are real parameters. Then we find a solution 

Cd = -11.1923, 

K = -0.01433. (35) 

This means we could reproduce the low-energy experimental data related to the charged 
fermion sector. As discussed in the previous section the Yukawa coupling matrices, Yio and 
Y1205 cire related to the corresponding mass matrices Mio and M120 such that 

Yio = ^ M,o, 
a"f sm p 

Yuo = . M120. (36) 

p^v sm p 



Here and /3" are the Higgs doublet mixing parameters introduced in ^13j, which are 
restricted in the range + < 1. Although these parameters are irrelevant to fit the 
low-energy experimental data of the fermion mass matrices, there are theoretical lower bound 
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on them in order for the resultant Yukawa coupling constant not to exceed the perturbative 
regime. Now we show one example of the Yukawa coupling matrices with fixed a" = 0.707, 

0.00304 0.00644 -0.00236 
Fio = I 0.00644 0.0235 -0.0502 | , (37) 
-0.00236 -0.0502 0.761 

O -0.000170 i -0.00557 i 
Yi2o=\ 0.000170 z O -0.0000226 i |. (38) 

0.00557 i 0.0000226 i O 

For the proton decay analysis, we search for the parameters which cancel the proton decay 
rate, completely. In the foUowings, we restrict the region of the parameters in the range 
(a")^ + (/5")^ = 1 (we assume a" and (3^ real for simplicity) . Then we actually find a 
solution that cancels the proton decay p — >■ K~^V. The resultant 2x2 color-triplet mass 
matrix ^13j has the foUowing form: 

„ „ / -0.0814 + 0.113i -0.872 -0.470i \ 
Mc-MGhx^ 0.872 - 0.4702 -0.0814 - 0.113^ J' 

For about Model 2, we have already found a solution in [13]. Therefore we just present the 
result here 

K = -0.00675 + 0.000309 i , 

\cd\ = 5.99 , (40) 

with tan /3 = 2.5. 



4 Neutrino oscillation data 

In our scheme adopted here, the right-handed neutrino mass matrix is left free from fitting 
the charged fermion data. We know the definite structure for the Dirac neutrino mass matrix. 
Therefore, by making use of the neutrino oscillation data, we can give a definite prediction 
for the right-handed neutrino mass matrix. 

Mr = MdM-^mI 

= MDUl^^dmg{m^\m^\ml^)U;^^^Ml. (41) 
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Here Uuns is the Maki-Nakagawa-Sakata (MNS) lepton mixing matrix, and in the Standard 
parametrization, it can be written as 



-iS 

C13C12 C13S12 SiaC 

^MNS = I — C23S12 — ■523Ci2Si3e"' C23C12 — S23Si2Si3e*'' S23C13 
S23S12 - C23Ci2Si3e*'^ -S23C12 - C23Si2Si3e*'' C23C13 

X diag(l,e^^e^^), (42) 

where Sij := sinOij, Cij := cos%, 6, /3, 7 are the Dirac phase and the Majorana phases, 
respectively. Now we take input the center values of a global analysis for the neutrino 
oscillation parameters after the recent KamLAND data [H], [fl 





= Am2i 


= 2.1 X 10~2 eV^ 


sin^ 6*® 


= 0.5, 




Am|, 


= ^"^ii 


= 8.3 X 10"^ eV^ 


sin^ 9q 


= 0.28, 




lUesl 


< 0.15. 





(43) 

Note that we can take both signs of Amgj^, Amg^ > O or Amg^ < 0. The former is called 
"normal hierarchy" (NH), the latter is called "inverted hierarchy" (IH), and we study both 
cases. Writing the lightest neutrino mass eigenvalue as m^, we can write the mass eigenvalues 
as 

mi = mi, 



m2 = \ rnj + Am| 



= Jmj + Aml + Aml, (44) 



for the case of NH, and 



mi = \/rnj + Am^ + Am|, 



m2 = ym^ + Am^, 

ms = mi, (45) 



for the case of IH. 



11 



Because we know the Dirac mass matrix M d, by inputting the above neutrino oscillation 
data, we can predict the mass matrix of the right-handed neutrinos as a function of Majorana 
phases, (3 and 7. As a result, we can have a prediction on the phenomena relating the right- 
handed neutrinos, such as the leptogenesis, lepton flavor violating processes (LFV), etc. 

Now we turn to the discussions about the baryon asymmetry of the universe based on 
the leptogenesis scenario [15]. In the leptogenesis scenario, lepton asymmetry is generated 
by the out-of-equilibrium decays of the right-handed neutrinos. The lepton asymmetry of 
the right-handed neutrino Ni is defined as 



W) 



At tree level, the decay width of Ni is given by 

{MdMI)u 



wy 



(46) 



r, 



M. 



Svrf 2 sin^ /3 

The CP asymmetry is produced for the first time on one-loop level as 
1 



(47) 



DM\))ii ^ 



Ml 



DJiji 



M2 



(Ml 



(4J 



where f{x) and gj{x) denote the one loop contribution from the vertex and the self-energy, 
respectively [T6] . 



In 



l + x 



X 



{x-l 



{x-iy+ {MdMI,),,/{8 vTf 2 sin (3) X 



■ (49) 



Once we are able to know the values of the lepton asymmetry e», we have to solve the 
Boltzmann equations in order to get the actual baryon asymmetry, rj = ns/n^- But we can 
use an approximate formula [T7|, 



3 X 10 



-3 



10-^ eV 



nii 



log 



mi 



10-3 eV 



(50) 



where 



mi 



{MnMl)n 
Ml 



(51) 



This formula is within a sufficiently good approximation for 10 ^ [e V] < mi < 10^ [eV]. 
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Now we we show the results of the baryon asymmetry of the universe in case of NH, 
in Fig. [H [2], and [31 The vertical axis represents the predicted baryon asymmetry of the 
universe 17 as a function of the hghtest neutrino mass mi [e V] in Fig. [T], as a function of |?7e3| 
in Fig. [2] and as a function of 6 [rad] in Fig. [3l In Fig. [H we have fixed \Ue3\ = 0.15 and the 
Majorana phases /3 and 7 as /5 = 7 = 7r/4. In Fig. [21 mi is fixed as mi = 10^^ [e V] and 
the Majorana phases are the same as Fig. [H For Fig. [H and [21 three lines (red, blue, green) 
represent the values 6 = n/A, tc/2, n, respectively. In Fig. [3l we have fixed mi = 10~^ [eV] 
and the Majorana phases (3 and 7 as /3 = 7 = 7r/4 but varied \Ue3\ = 0.15, 0.10, 0.05, 
which correspond to the three lines (red, blue, green). These results show that Model 1 can 
not reproduce the observed baryon asymmetry in the Standard leptogenesis scenario, though 
all possible values of free parameters have not been exhausted. But in strictly speaking, 
any (local) supersymmmetric models have a serious problem, so called, gravitino problem 
[T8l [T9l [20] . That says if we take a naively expected value for the garavitino mass 1713/2 of 
order the weak sacle m3/2 ~ 100 [GeV] the lifetime is shorter than 1 [ SGc] , cis cL result the 
primordial gravitino dacays after the big-bang nucleosynthesis (BBN). Recently, the updated 
analysis for the hadronic decay of the gravitino shows that the reheating temperature is very 
strictly constrained as Tr < 10^~^ [GeV] [20]. Therefore the results of our Model 1 for the 
baryon asymmetry of the universe should not taken so seriously. 

On the other hand, in Fig. [11 [5l [6] and [7] we show the results of the baryon asymmetry of 
the universe in case of Model 2. The vertical axis represents the predicted baryon asymmetry 
of the universe 77 as a function of the CP-phase a [rad]. In Fig. [H we have fixed mi = 10"'^ 
[eV] and |?7e3| = 0.15, and Majorana phases /5 and 7 as /5 = 7 = 7r/4. In Fig. [5l we have 
adopted the same parameters as Fig. [l]except for the Majorana phases /? = 7 = n/2. 
Fig. [6] is the zoomed up picture of the experimentally allowed parameters region of Fig. 
[S] for the Dirac CP-phase S = n /4. Fig. [7] is the same diagram as Fig. [H except for the 
Majorana phases /3 = 7 = vr. You can see that this model can reproduce the observed 
baryon asymmetry within the Standard leptogenesis scenario. 

Since Model 2 has passed the test for the leptogenesis, it is possible to go further steps, 
e. g. analysis of the proton lifetime r(p — K~^i?) in a mode p —>■ K^v. The predicted 
values as a function of a [rad] are depicted in Fig. [S] Though in general we can vary 
the parameters in the Higgs potential as possible, we pick up one solution satisfying the 
experimental constraint. In concrete, we used the 2x2 color-triplet Higgsino mass matrix 
[13] of the following form. 



which corresponds to the parameters with 9 = (p = ip' = 7r/4in the notation of 
Ref. [13]. The horizontal line indicates the current experimental lower bound. Therefore, 




(52) 
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these results show that there exists a parameter region which simultaneously reproduces a 
required baryon asymmetry of the universe and satisfies the current experimental bound on 
the proton lifetime: 



5 Conclusion 

We have proposed a phenomenologically viable SO(IO) grand unification model which use a 
set of Higgs superfields 10 + 16 + 16 + 45. In this framework, the data fitting of the charged 
fermion mass matrices have been re-examined. This model indeed can fit the low-energy 
experimental data relating the charged fermion masses and mixings. As for the neutrino 
sector, we have inputted the neutrino oscillation data to constrain right-handed Majorana 
neutrino mass matrix. Then the unknown parameters are used to fit the leptogenesis and 
for the proton decay. In such a framework, two models, Model 1 and Model 2, have been 
considered. The former is effectively equal to the model with the Yukawa couplings to the 
10 + 120 Higgses and the latter to the model with the Yukawa couplings to the 10 + 126 
Higgses with respect to charged Fermions. It should be remarked in Model 2 that has 
the different origin from that of the charged fermion mass matrices unlike the conventional 
minimal SO(IO) model [6] [7j. Then it has been found that our model (Model 2) is consistent 
with all experimental constraints, especially, the required baryon asymmetry of the universe 
?7 can be produced, keeping in mind the proton lifetime constraint. Though Model 1 does 
not produce the sufficient baryon asymmetry of the universe within the Standard framework 
of the leptogenesis, there is a serious problem in any SUSY models, so called, gravitino 
problem [TU [191 120] • The decay products of the primodial gravitino destroy the produced 
baryon asymmetry of the universe and need further study to consider some extensions of the 
standad leptogensis scenario, such as Affleck-Dine scenario [2T] . 
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Figure 1: The baryon asymmetry of the universe 77 as a function of the lightest neutrino mass 
mi [e V] for Model 1. Three lines (red, blue, green) correspond to the values 5 = 7r/4, 7r/2, vr, 
respectively. In this figure we have fixed jt/esl = 0.15 and the Majorana phases (3 and 7 as 
/5 = 7 = 7r/4. 
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Figure 2: The baryon asymmetry of the universe rj as a function of {Uesl Model 1. Three 
lines (red, blue, green) correspond to the values 5 = 7r/4, tt/2, tt, respectively. In this figure, 
we have fixed mi = 10~^ [e V] and the Majorana phases (3 and 7 as /5 = 7 = n /4. 
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Figure 3: The baryon asymmetry of the universe 77 as a function of 6 [rad] for Model 1. Three 
lines (red, blue, green) correspond to the values jf/esl = 0.15, 0.10, 0.05, respectively. In this 
figure, we have fixed mi = 10~^ [e V] and the Majorana phases (3 and 7 as /5 = 7 = n /4. 
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Figure 4: The baryon asymmetry of the universe 77 as a function of a [rad] in Model 2. 
Three lines (red, blue, green) correspond to the values 5 = 7r/4, 7r/2, tt, respectively. In 
this figure, we have fixed mi = 10"^ [eV] and iL^esI = 0.15, and the Majorana phases as 
/5 = 7 = 7r/4. 
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Figure 5: Same as Fig. HJ but for the Majorana phases /3 



7 



7r/2. 



6. 


9 


•10 


-10 


6. 


8 


•10 


-10 


6. 


7 


•10 


-10 


^ 6. 


6 


•10 


-10 


6. 


5 


•10 


-10 


6. 


4 


•10 


-10 


6. 


3 


•10 


-10 



o 



2 3 4 

o [ rad] 



5 



Figure 6: The experimentally allowed region in Fig. [5]are zoomed in for 5 = n 
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Figure 7: Same as Fig. HJ but for the Majorana phases /5 = 7 = tt. 
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Figure 8: The proton lifetime r(p —>■ K^v) in a mode p — *• K^v as a function of a [rad] for 
Model 2. Though in general we can vary the parameters in the Higgs potential as possible, 
we pick up one solution satisfying the experimental constraint. The horizontal line indicates 
the current experimental lower bound. 
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